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ABSTRACT 

O 

^ \ Using three-dimensional resistive MHD nested grid simulations, we investigate 

the driving mechanism of outflows and jets in star formation process. Starting 
with a Bonnor-Ebert isothermal cloud rotating in a uniform magnetic field, we 
calculated cloud evolution from the molecular cloud core (n c = 10 4 cm -3 , r = 
4.6 x 10 4 AU) to the stellar core (n c = 10 22 cm 3 , r ~ 1-R©), where n c and r denote 
the central density and radius of each object, respectively. In the collapsing cloud 
core, we found two distinct flows: Low- velocity outflows (~5kms _1 ) with a wide 
opening angle, driven from the adiabatic core, and high- velocity jets (~30 km s" 1 ) 
with good collimation, driven from the protostar. High-velocity jets are enclosed 
by low-velocity outflow. The difference in the degree of collimation between the 
two flows is caused by the strength of the magnetic field and configuration of the 
magnetic field lines. The magnetic field around an adiabatic core is strong and has 
an hourglass configuration; therefore, the low- velocity outflow from the adiabatic 
core are driven mainly by the magnetocentrifugal mechanism and guided by the 
hourglass-like field lines. In contrast, the magnetic field around the protostar is 
weak and has a straight configuration owing to Ohmic dissipation in the high- 
density gas region. Therefore, high-velocity jet from the protostar are driven 
mainly by the magnetic pressure gradient force and guided by straight field lines. 
Differing depth of the gravitational potential between the adiabatic core and the 
protostar cause the difference of the flow speed. Low- velocity outflows correspond 
to the observed molecular outflows, while high-velocity jets correspond to the 
observed optical jets. We suggest that the protostellar outflow and the jet are 
driven by different cores, rather than that the outflow being entrained by the jet. 



1. Introduction 

The observations indicate that outflow is ubiquitous in the star formation process, and 
flows from the protostars have two or more distinct velocity components (El; Q). Typically, 
a flow from the protostar is composed of a low- velocity component (LVC) with 10-50 km s _1 
and a high- velocity component (HVC) with ~100kms _1 . There is a clear trend toward 
higher collimation at higher flow velocity (|3|). Since flows from the protostar have various 
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morphological and kinematical properties, they cannot be explained by a single-class model. 
The flows that originated from the protostar are typically classified into two types: molecular 
outflow observed mainly with CO molecules (0), and optical jet observed by optical emission 
(Q). Molecular outflows observed by CO line emission exhibit a wide opening angle (j^l) and 
slower velocity of 10 — 50kms _1 (J3|), while optical jets observed by optical emission exhibit 
good collimation and higher velocity of 100 — 500 km s -1 (jsl). Observations indicate that 
around each protostar, high-speed jets with a narrow opening angle are enclosed by a low- 
velocity outflow with a wide opening angle (0). However, the driving mechanism of these 
flows are still unknown. 

In this study, we calculate cloud evolution from the molecular cloud core (n c = 10 4 cm -3 , 
r c = 4.6 x 10 4 AU) to stellar core formation (n c ~ 10 22 cm -3 , r c ~ 1-R©) using three- 
dimensional resistive MHD nested grid method, study the formation process of jets and 
outflows, and show the driving mechanisms of these flows. 



Model 



Our initial settings are almost the same as those of (@; BE Q 0)- We solve the 



resistive MHD equations including self-gravity (see, eqs.1-5 of (lid )). We adopt a spherical 
cloud with critical Bonnor-Ebert density profile having pbe = 3.841 x 10 _20 g cm -3 (n^ = 
10 4 cm -3 ) of the central (number) density as the initial condition. The critical radius for a 
Bonnor-Ebert sphere R c = 6.45 c s [47tG/9be,o] corresponds to R c = 4.58 x 10 4 AU for our 
settings. Initially, the cloud rotates rigidly (Qq = 7 x 10 15 s _1 ) around the z-axis and has 
a uniform magnetic field (B = 17/iG) parallel to the z-axis (or rotation axis). To promote 
contraction, we increase the density by 70% from the critical Bonnor-Ebert sphere. The 
initial central density is therefore po = 6.53 x 10~ 20 gcm -3 (no = 1.7 x 10~ 4 cm~ 3 ). 



We adopt the nested grid method ()8|; [13|; ll4j) to obtain high spatial resolution near 
the center. Each level of a rectangular grid has the same number of cells (64 x 64 x 32), 
although the cell width h(l) depends on the grid level I. The highest level of a grid changes 
dynamically. The box size of the initial finest grid I = 1 is chosen to be 2 4 R C , where 
R c denotes the radius of the critical Bonnor-Ebert sphere. A new finer grid is generated 
whenever the minimum local Jeans length Aj becomes smaller than 8 h(l max ). The maximum 
level of grids is restricted to / max = 30. 
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Fig. 1. — Bird's-eye view. The structure of high-density region (p > 0.1p c ; red iso-density 
surface), and magnetic field lines (black-and-white streamlines) are plotted in each panel. 
The structures of the jet (v > 7kms _1 ) and outflow (v > 0.5 km s -1 ) are shown by iso- 
velocity surfaces, respectively. The density contours (false color and contour lines), velocity 
vectors (thin arrows) on the mid-plane of x — 0, y — 0, and z =0 are, respectively, projected 
in each wall surface. 
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3. Results 

The molecular gas obeys the isothermal equation of state with temperature of ~ 10 
K until n c ~ 5 x 10 10 cm -3 (isothermal phase), then cloud collapses almost adiabatically 
(5 x 10 10 cm -3 < n c < 10 16 cm -3 ; adiabatic phase) and quasi-static core (i.e., first core) 



forms during the adiabatic phase (Il5l ; Il6l ). In our calculations, the first core forms when 



the central density reaches n c ~8x 10 cm . The magnetic flux is removed from the first 



core during the adiabatic phase by the Ohmic dissipation (1171 ) . After central density reaches 
n c ~ 10 16 cm" 3 , the equation of state becomes soft reflecting the dissociation of hydrogen 
molecules at T ~ 2 x 10 3 K, and collapses rapidly. By this epoch, the central temperature 
becomes so high that the thermal ionization of Alkali metals reduces the resistivity and so 
that Ohmic dissipation becomes ineffective. Thus, the magnetic field becomes strong again 



as central region collapses. The second core (or protostar) (I la ) forms at n c c± 10 21 cm -3 . 
The magnetic field strength increases rapidly after the second core formation epoch (n > 
10 21 cm -3 ), because the shearing motion between the second core and ambient medium 
amplifies the toroidal magnetic field around the second core. 

Figure [1] shows the structure of the low- and high- velocity flows, and the configuration 
of the magnetic field lines. It also shows the shapes of the first core (left panel; the projected 
density contours on the wall) and the protostar (right panel; the red isosurface). The purple 
and blue surfaces in Figure [T] indicate the iso- velocity surface of v z = 5 km s" 1 and v z = 
0.5 km s -1 , respectively. The flow inside the purple iso-velocity surface has a velocity of 
v > 0.5 km s -1 (low- velocity component; LVC), while the flow inside the blue iso-velocity 
surface has a velocity of v > 5kms _1 (high- velocity component; HVC). The HVC is enclosed 
by the LVC. The LVC flow is mainly driven from the first core, while the HVC flow is mainly 
driven from the protostar. The LVC and HVC are strongly coiled by the magnetic field lines 
anchored to the first core and protostar, respectively. 

Our results show that the flow appearing around the first core has a wide opening angle 
and slow speed, while the flow appearing around the protostar has a well-collimated structure 
and high speed, as shown in Figure [2J The speed difference is caused by the difference of 
the depth in the gravitational potential. The flow speed corresponds to the Kepler speed of 
each object. Because the first core has a shallow gravitational potential, its flow is slower. 
The flow driven from the protostar, which has a deeper gravitational potential, has a high 
speed. In our calculations, the low- and high- velocity flows have speeds of i>lvf — 3kms _1 
and -uhvf — 30kms _1 , respectively. These speeds are slower than those of observations. 
Typically, observed molecular outflow and optical jet have speeds of f ou t,obs — 30kms _1 , and 
"Wjet.obs — 100 km s -1 , respectively. However, since the first and second cores (protostar) have 
mass of Mf c = 0.01 M Q and M sc ~ 1O~ 3 M , respectively, at the end of the calculations, 
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each core increases its mass in the gas accretion phase. The Kepler speed increases with the 
square root of the mass. When the mass of each core increases by 100 times, the Kepler 
speed increases 10 times. Thus, the speed of the low- and high-velocity flows may increase by 
10 times, and reach i> L vf — SOkms^ 1 and wrvf = 300 km s -1 , respectively, which correspond 
to typical observed values. 

4. Discussion 

Observation shows that the molecular outflows have wide opening angles and low flow 
speeds, while the optical jets have good collimation and high flow speeds. Molecular outflow 
has been considered to be entrained by the optical jet driven from a circumstellar disk around 
the protostar. In this study, we calculated the cloud evolution from the molecular cloud core 
to protostar formation, and found that two distinct flows are driven from different objects, 
and the observed features of molecular outflow and optical jet were naturally reproduced. 
Thus, we expect that the low-velocity flow from the first core corresponds to the molecular 
outflow, while the high-velocity flow from the protostar corresponds to the optical jet. 

The different collimation of low- and high-velocity flow is caused both by the config- 
uration of the magnetic field lines around the drivers and their driving mechanisms. The 
magnetic field lines around the first core have an hourglass configuration because they con- 
verge to the cloud center as the cloud collapses, and Ohmic dissipation is ineffective before 
the first core formation. In addition, the centrifugal force is more dominant than the Lorentz 
force in the low- velocity flow (molecular outflow). Thus, the flow appearing near the first 
core is mainly driven by the magnetocentrifugal wind mechanism (disk wind). On the other 
hand, near the protostar, the magnetic field lines are a straight, and the magnetic pressure 
gradient mechanism is more effective for driving the high- velocity flow (optical jet). The 
magnetic field lines straighten by the magnetic tension force near the protostar because 
the magnetic field is decoupled from the neutral gas. However, the magnetic field lines are 
strongly twisted in the region in close proximity to the protostar, where the magnetic field 
is coupled with the neutral gas again. Thus, the strong toroidal field generated around 
the protostar can drive the high- velocity flow (optical jet), which is guided by the straight 
configuration of the magnetic field. 

Our calculations do not completely reject the well-known concept that the observed 
molecular outflow is entrained by the optical jet, because we calculate the formation of the 
jet and outflow only in the early star-formation phase. Further long-term calculations are 
needed to understand the mechanism of the optical jet and molecular outflow in more detail. 
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jet (high velocity) 




Fig. 2. — Schematic view of the jet and outflow driven from the protostar and the first core, 
respectively. 
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